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Abstract 

We have carried out the basic research for the accelera- 
tor and tokamak control system based on the Experimental 
Physics and Industrial Control System (EPICS). We have 
used the process database and the state notation language 
(SNL) in the EPICS to develop the simulator which repre- 
sents as a virtual machine. In this paper. In this paper, we 
introduce the simulator of the global orbit feedback system 
as an example. This simulates the global orbit feedback 
system under the constraint conditions for Pohang Light 
Source (PLS) storage ring. We describe the details of the 
feedback algorithm and the realization of the simulator. 

1 INTRODUCTION 

When the beam position change takes place during the 
conventional global orbit correction processes, the photon 
beam through the beamline is affected, and it results in the 
alignment of mirrors and monochromators. This is par- 
ticularly severe for a long beamline such as the undulator 
beamlines. This problem can be overcome by introduc- 
ing a local bump at the particular beamline. However, for 
some light sources, there are not enough corrector magnets 
to generate local bumps as much as needed. This difficulty 
can be overcome when we correct the COD under the con- 
dition where the beam positions at particular points are not 
changed. This is our main objective to develop a method of 
the closed orbit correction under constraint conditions. 



2 THEORY 

2.1 Ordinary COD correction and regulariza- 
tion 

The beam position is normally described as a vector |x) 
measured by M beam position monitors (BPM). In order to 
correct the COD, we need N corrector magnets with then- 
strengths described as a vector |k). When the corrector 
magnets kick a beam, the new beam positions |y) can be 
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described as follows. 



|y)=R|k) + |x). 



(1) 



Here, R is called the response matrix of (M x N) dimen- 
sions whose components are given by 



R, 



2 sin m> 



COS (|VP, - - TTV) , 



(2) 



where v is the betatron tune of the storage ring, and (/% , \l/ j) 
and (f3j, tyj) are the beta function and the phase function 
for the i th BPM and j th corrector magnet, respectively. In 
order to reduce the COD, we have to choose the kick of 
each corrector magnet satisfying 



R T R |k) + R T |x) = 



(3) 



It is called the PSINOM algorithmic . 

The COD correction is actually a minimization proce- 
dure of S defined as 



S 



{ (k| R T R |k) + 2 (x| R |k> + (x |x) } 



(4) 



By using the relations of vector operators which are hyper- 
dimensional gradient operators, we can get the same result 
of PSINOM algorithm as shown in eq. (fy. 

This algorithm includes an inversion procedure of the 
matrix R T R. In some cases, we can get unacceptable cor- 
rections due to the ill-posedness of R T R. A regularization 
method is introduced to avoid this problem. In this case, S 
is written as 



S = i { (k| R T R |k) + 2 (x| R |k) 

1 



f (x |x}} 
<k|a|k). 



(5) 



where a is the regularization parameter. Then the mini- 
mum corrector kicks can be determined by 



(R T R + al) |k)+R T |x} = 



(6) 



This equation represents the modified PSINOM algorithm, 
and we can relax the inversion problem of the singular 
matrix by using the diagonal matrix al when R T R is 
singular [0]. 



2.2 Method with constraint conditions 

After having the new closed orbit with the minimum distor- 
tion from eq. (^), the new orbit is generally different from 
the original orbit. Sometimes, this difference can be taken 
place at very sensitive locations such as the entrance and 
the exit of an undulator. If the beamline is well aligned for 
this undulator, a COD correction should be avoided in this 
region. A constraint condition can be described in terms of 
the beam position at i th BPM such as 

(Ri |k) + x 0i = Xi . (7) 

Here, (R; | is the i th row of the response matrix R. Also, 
XQi and Xi are the beam positions before and after the cor- 
rection, respectively. Since we want to keep this position 
unchanged, (R, |k) should be zero. If there are L BPMs 
involved in the constraint condition, we can write the con- 
straint condition as follows. 

C T |k) = 0. (8) 

Here, C T is the (L x N) sub-matrix of the response matrix. 
Each component of C T corresponds to the BPM involved 
in the constraint condition. We also assume that |k) has a 
non-trivial solution. 

We now add this constraint condition to the modified 
PSINOM algorithm to obtain the new S such as, 

S = i [(k| R T R |k> + 2 (x| R |k) + (x|x>] 

+~ <k| a |k> + <r| C T |k) . (9) 

Here, (r| is the Lagrangian multiplier, and it is an L dimen- 
sional vector. By following the derivative to the corrector 
strength, we can get the vector |k) which minimizes the 
closed orbit distortion outside the constraint region such 

as, 

(A + aI)|k)+R T |x) + C|r) =0. (10) 

Here, we define the square matrices of N x N dimensional 
A and L x L dimensional D as follows. 

A = R T R, (11) 
D = C T (A + cd) -1 C . (12) 

Eq. (|Io|) can be rewritten as 

|r> = -D 1 ^ (A + air 1 |x> • (13) 

Now, we can remove the Lagrangian multiplier \T) in eq. 
(|I(J) by using the above equation. Then, we can finally get 
the kick values of the corrector magnets as follows. 

|k) = -(A + al) -1 {i-CD^C 7 (A + al) -1 } 

R T |x) . (14) 



3 DEVELOPMENT OF SIMULATOR 
USING EPICS 

The algorithm developed in the previous section has been 
successfully tested for the Pohang Light Source (PLS) op- 
eration. The new correction code is written in C language 
and it is installed in one of the operator consoles which 
is a SUN workstation. Although the PLS control system 
is not using Experimental Physics and Industrial Control 
System (EPICS) which is used in many accelerator labo- 
ratories world-widely, there is a plan to upgrade the con- 
trol system based on EPICS technology^ Q]. As a part of 
such upgrade activity, we have started the development or- 
bit correction algorithm with EPICS. Since we do not have 
any EPICS-based control system yet, we decide to develop 
the orbit correction simulator using EPICS. 

In order to seek the way to adopt our orbit correction al- 
gorithm into EPICS, we have considered two ways: one is 
based on the subroutine record and the other is using the 
state notation language (SNL) program. The first method 
needs new record support that runs the orbit correction al- 
gorithm. To do this, the correction code is required to be 
written upon the protocol required by the record support 
such as the entry structure and the callback structure. This 
approach gives relatively fast response because this method 
uses database access and can access the record by process 
passive mode[Q, ^]. This is a good feature for the realtime 
system but the large portion of the orbit correction code 
we have already tested must be rewritten according to the 
protocol of the record support. 

One of important tasks of EPICS input/output controller 
(IOC) is the sequencer that runs programs written in SNL. 
The SNL considers the control object as the state machine 
and treats transitions between states. The sequencer mon- 
itors the transitions for the SNL and runs callback func- 
tions using the entry table the corresponding program writ- 
ten in SNL. Since the sequencer accesses to the record via 
channel access (CA) and the record access is only possible 
through non-process passive mode, there are some restric- 
tions in access time or the treatment of records. However, 
this method can directly imbed the program written in C 
language. Also, unlike the subroutine record, this method 
can remove program tasks without rebooting the system, 
which gives the code debugging very easy[Q, g], Thus, the 
second method gives more benefits when the system does 
not require heavy realtime demands. Upon reviewing the 
two methods, we have decided to use the latter method. 

3. 1 SNL program 

The orbit correction simulator is developed in IOC level 
and has the SNL program imbedding C codes and several 
database records, as shown in Fig. [y. There are two parts 
in the SNL program: one for the feedback including the 
orbit correction algorithm and the other for the simulator 
which emulates the PLS storage ring. The latter calculates 
orbit changes from ai:KICK$(kickjio) records which store 
the corrector strengths obtained from the response matrix 




Figure 1 : Schemematic diagram of the orbit correction simulator 



measured at the PLS. This process is actually a product of 
matrix and column vectors. The calculated orbit changes 
are stored at ai:CorrOrbit$(bpmjio) records. In the next 
step, the corrector strengths s calculated by the algorithm 
described in the previous section and the results are stored 
atal:KICK$(kick.no). 

3.2 Database and record linkage 

The SNL program we have developed consists of four state 
sets. Since they are linked together and interacted dy- 
namically via EPICS database, we need to understand the 
database and its records as well as the linkage between 
them. The ai:KICK$(kickjio) record uses the analog in- 
put record to represent the corrector strength. There are 
70 records altogether as representing 70 correctors in the 
PLS storage ring. They are linking the simulator and the 
feedback in a parallel manner. The index $(kickjio) is the 
integer value between one and 70. 

The records of ai:OrgOrbit$(bpmjio) and 
ai:CorrOrbit$(bpmju>) are the orbit distortion before 
the feedback is applied and the orbit change after the 
feedback, respectively. They are using analog input 
records. Since these records represents the changes at 
the PLS BPM, $(bpmjio) is the interger value between 
one and 108. The superposition of these two records 
gives new orbit. This can be done by the calculation 
record calc:Orbit$(bpmjio). Since this calculation 
record is linked forwardly from ai:OrgOrbit$(bpmjio) 
and ai:CorrOrbit$(bpmjio), the superposition is newly 
calculated whenever the values of two records are changed. 

On the other hand, aimumBPM, aimumCorrector, 
ai. startCONSTR, and ai. endCONSTR are representing the 
BPMs and correctors used in the correction algorithm, start 
and end point of the constraint region, respectively. They 
are using analog input records and notifying the SNL pro- 
gram if necessary. 

The other records are binary input records and they are 



representing necessary status flags. They are used as the 
mediator of state transitions between state sets in the SNL 
program. 

4 CONCLUSION 

We have developed the COD correction algorithm under 
the constraint condition where the beam position at partic- 
ular point is not changed. The new algorithm is based on 
the modified PSINOM algorithm which includes the regu- 
larization process in order to avoid the inversion problem 
of the ill-posed response matrices. 

We have confirmed that this algorithm is working well 
and is in good agreement with the experimental results [Q]. 
Even though the PLS is planning to upgrade their control 
system with EPICS, there is no working EPICS based con- 
trol system at PLS. Due to this, we have developed the orbit 
correction simulator using C-code embedded SNL program 
based on EPICS technology. This simulator part can be re- 
placed by the real control system with minor changes after 
the completion of the upgrade. 
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